Rotating skyrmion lattices by spin torques and field or temperature gradients 
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Chiral magnets like MnSi lorm lattices of skyrmions, i.e. magnetic whirls, which react sensitively 
to small electric currents j above a critical current density jc- The interplay of these currents with 
tiny gradients of either the magnetic field or the temperature can induce a rotation of the magnetic 
pattern for j > jc- Either a rotation by a finite angle of up to 15° or - for larger gradients - 
a continuous rotation with a finite angular velocity is induced. We use Landau-Lifshitz-Gilbert 
equations extended by extra damping terms in combination with a phenomenological treatment 
of pinning forces to develop a theory of the relevant rotational torques. Experimental neutron 
scattering data on the angular distribution of skyrmion lattices suggests that continuously rotating 
domains are easy to obtain in the presence of remarkably small currents and temperature gradients. 

PACS numbers: 



INTRODUCTION: SPINTORQUES AND 
SKYRMION LATTICES 



Manipulating magnetic structures by electric current 
is one of the main topics in the field of spintronics. By 
strong current pulses one can, for example, switch mag- 
netic domains in multilayer deviceaii^, induce microwave 
oscillations in nanomagnets^ or move ferromagnetic do- 
main walls^. The latter effect may be used to develop 
new types of non- volatile memory devices^. It is there- 
fore a question of high interest to study the coupling 
mechanisms of currents to magnetic structuresii^. 

Here, the recent discoverjiSiiS of the so-called skyrmion 
lattice in chiral magnets like MnSi provides a new oppor- 
tunity for studying the manipulation of magnetism by 
electric currents both experimentally and theoretically. 
The skyrmions in MnSi form a lattice of magnetic whirls, 
similar to the superfiuid whirls forming the vortex lattice 
in type-II superconductors. While in ordinary ferromag- 
nets, currents couple only to the canted spin configura- 
tions at domain walls, the peculiar magnetic structure 
of the skyrmion lattice allows for an efficient hulk cou- 
pling. Furthermore, the smooth magnetic structure of 
the skyrmion lattice decouples efficiently from the un- 
derlying atomic lattice and from impurities. As a conse- 
quence, it was observedi^ that the critical current density 
needed to affect the magnetic structure was more than 
five orders of magnitude smaller than in typical spin- 
torque experiments. 

These low current densities open opportunities for new 
types of experiments to study quantitatively the physics 
of spin transfer torques. Due to the much lower cur- 
rent densities it is now possible to perform spintorque 
experiments in bulk materials and thus avoid the surface 
effects that dominate in nanoscopic samples. Moreover 
for smaller currents the effects of heating and Oersted 
magnetic fields created by the current are suppressed. 
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Figure 1: Schematic plot of the forces on a skyrmion lattice 
perpendicular and parallel to the current fiowing in vertical 
direction. For a static, non-moving skyrmion lattice the red 
horizontal arrows correspond to the Magnus force and the 
green vertical arrows to dissipative forces. In the presence of 
a temperature or field gradient, these forces change smoothly 
across a domain, thereby inducing rotational torques which 
depend sensitively on the relative orientation of current and 
gradient (and on the direction in which the skyrmion lattice 
moves). Small black arrows: local orientation of the magneti- 
zation projected into the plane perpendicular to the magnetic 
field B. In each unit cell the magnetization winds once around 
the unit sphere. 



In this paper we suggest experiments and develop a 
theory with the goal to exploit the rotational motion in- 
stead of just translational motion to investigate the inter- 
play of electric currents and moving magnetic structures. 
Our theory is directly motivated by recent experiments^^, 
where a change of orientation of the skyrmion lattice as a 
function of the applied electric current was observed with 
neutron scattering. In Ref. (lo| we have shown that the 
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rotation arises from the interplay of a tiny thermal gradi- 
ent parallel to the current and the Magnus forces arising 
from the spintorque coupling of current and skyrmion lat- 
tice. For example, the rotation angle could be reversed 
by reversing either the current direction or the direction 
of the thermal gradient. 

The basic idea underlying the theoretical analysis of 
our paper is sketched in Fig. [TJ In the presence of an 
electric current several forces act on the skyrmion lat- 
tice. First, dissipative forces try to drag the skyrmion 
lattice parallel to the (spin-) current. Second, the inter- 
play of dissipationless spin-currents circulating around 
each skyrmion and the spin-currents induced by the elec- 
tric current lead to a Magnus force oriented perpendic- 
ular to the current for a static skyrmion lattice (for the 
realistic case of moving skyrmions the situation is more 
complicated). In the presence of any gradient across the 
system (e.g. a temperature or field gradient), indicated 
by the color gradient, these forces will vary in strength 
across a skyrmion domain. 

As in the experiment, we assume that the gradients 
are tiny: on the length scale set by the skyrmion distance 
the gradients have negligible effects. However, multiply- 
ing the tiny gradient with a large length, i.e., the size of 
a domain of the skyrmion lattice (which can be^i several 
hundred nm), one obtains a sizable variation of the forces 
across the domain. These inhomogeneous forces can give 
rise to rotational torques. Whether the torque arises from 
the Magnus forces or the dissipative forces depends, how- 
ever, on the relative orientation of current and gradient 
and also on the direction in which the skyrmion lattice 
drifts. Fig.[T]gives a simple example: if, for example, cur- 
rent and gradient are parallel to each other (right panel) 
the forces perpendicular to the current direction (red hor- 
izontal arrows) give rise to rotational torques while the 
parallel forces do not contribute. The situation is re- 
versed when current and gradient are perpendicular (left 
panel). 

We therefore suggest to use the rotation of magnetic 
structures as a function of the relative orientation of cur- 
rent and further gradients as a tool to explore the cou- 
pling of magnetism and currents. We will show that the 
resulting rotations depend very sensitively both on the 
relative size of the various forces affecting the skyrmion 
dynamics and on how these forces depend on the induced 
gradients. While we apply our theory here to skyrmion 
lattices, our theoretical approaches can also be used for 
other complex magnetic textures and our results should 
also have ramifications for other setup a^^i^'^ . Quantita- 
tively, we will only study the role of gradients induced 
by changes in temperature or magnetic field but other 
options are also possible. For example, macroscopic vari- 
ations of the cross section of a sample will lead to gradi- 
ents in the current density. Also changes in the chemical 
composition or strain in the sample can induce gradients. 

It is also essential to investigate the effect of pin- 
ning of the magnetic structure by inhomogeneities arising 
from crystalline imperfections. Inhomogeneities distort 



the perfect skyrmion lattice and lead to forces prohibit- 
ing (up to a very small creep) the motion of the mag- 
netic structure as long as the current is below a criti- 
cal value, J < jc- Also for j > jc, inhomogeneities in- 
duce an effective, velocity dependent frictional force on 
the moving skyrmion lattice connected to local, time- 
dependent distortions of the skyrmion lattice. Pinning 
has widely been studied both experimentally and the- 
oretically for charge density waves and vortex lattices 
in superconductorsi^^— . As the dynamics of skyrmions 
differs qualitatively (and quantitatively) from these two 
cases it is not clear which of these results can be trans- 
ferred to skyrmion lattices. Due to the non-linear depen- 
dence of the pinning forces on the velocity, they can not 
be described by a simple damping term. Within this pa- 
per we will not try to develop a theory of pinning but will 
instead use a simple phenomenological ansatz to describe 
and discuss pinning effects. 

Rotational torques can also arise in the absence of the 
types of gradients discussed above. In Ref. [l3 we have 
studied the role of distortions of the skyrmion lattice by 
the underlying atomic lattice extending the methods used 
by Thielei^ to rotational torques (this method will also be 
used below). Such distortions indeed induce small rota- 
tional torques in a macroscopically homogeneous system, 
i.e. without any external gradients. Similarly, also distor- 
tions induced by disorder can induce rotational torques 
without external gradients as has been discussed in the 
seminal paper by Hanger and SchmicU^. But all these 
effects are very small and have not been observed in the 
experimental setup of Ref. [lo| as no rotation has been 
observed in the absence of gradients. Therefore they will 
be neglected in the following. 

In the following we will first describe briefly the rele- 
vant Ginzburg-Landau model and the Landau-Lifshitz- 
Gilbert equation used to model the dynamics of the 
skyrmions. Here we include a novel damping term a' 
recently introduced in Refs. [20||2T| (we also add the cor- 
responding /3' term). We then derive effective equations 
for the translational and rotational mode where pinning 
physics is taken into account by an extra phenomeno- 
logical term. This allows to develop predictions both for 
static rotations by a finite angle and continuous rotations. 
In the light of our results we interpret experimental re- 
sults on the angular distribution of skyrmion lattices in 
the presence of currents and gradients. 



II. SETUP 

A. Ginzburg-Landau model 

The starting point of our analysis is the standard 
Ginzburg-Landau model of a chiral magnet in the pres- 
ence of a Dzyaloshinskii-Moriya intcraction22i2^. After a 
rcscaling of the length r, the local magnetization M{r) 
and the magnetic field B the free energy functional re- 
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duces to^ 

F = y d'V [(1 + t)M^ + {VMf 

+2M- (V X M) + -S-M] , (1) 

Here t cx T—T^^ parametrizes the distance to the mean- 
field phase transition at B = from a phase with hehcal 
magnetic order [t < 0) to a paramagnetic phase {t > 
0) ^^!^'^ . In the presence of weak disorder t (and strictly 
speaking also the prefactors of all other terms) fluctuates 
slightly as a function of r. 

The skyrmion lattice (stabilized by thermal fluctua- 
tions) exists for a small temperature and field ranged. 
It is translationally invariant parallel to B and shows a 
characteristic winding of the magnetization in the plane 
perpendicular to B, see Fig. [TJ 

B. Landau-Lifshitz-Gilbert equation 

To describe the dynamics of the orientation fl{r,t) = 
M{r,t)/\M{r,t)\ of the magnetization M{r,t) in the 
presence of spin-transfer torques due to electric cur- 
rents we use the standard Landau-Lifshitz-Gilbert (LLG) 
equation,— iSii^ extended by a new dissipative term^SiSi 

n-(d^nx (dt + ^vsV)n)]d^n. (2) 

a' i 

Here Vg is an effective spin velocity parallel to the spin 
current density. More precisely, for smooth magnetic 
structures with constant amplitude of the magnetization 
it is given by the ratio of the spin current^ and the size 
of the local magnetization, \M\. In a good metal (for ex- 
ample, MnSi) Vs is expected to be parallel to the applied 
electric current and to depend only weakly on temper- 
ature and field. The magnetization precesses in the ef- 
fective magnetic field Weft ~ ^JT Strictly speaking 
Eq. ([2]) is only valid for a constant amplitude of the mag- 
netization, \M\ = const. Since \M\ varies only weakly^ 
in the skyrmion phase, we use as a further approxima- 
tion Hcs ~ ~JTJM where M is the average local 
magnetization, = (Af^). 

The last two terms in Eq. ([2]) describe dissipation. 
a is called the Gilbert damping and /3 parametrizes 
the dissipative spin transfer torque. The new damping 
term prop ortional to a' was introduced (for = 0) in 
Refs. \MM- It arises from the ohmic damping of elec- 
trons coupled by Berry phases to the spin texture as can 
be seen by rewriting Eq. ^ in the form 

- = Ain X (dt +VsV)n + aM (dt + -v^v) Cl 
sn \ a J 



where = Cl ■ {difl x dti^) can be interpreted as the 
emergent electric field and B^ = ^ajk^^ ■ {djft x 9fcf2) as 
the emergent magnetic fiel d^^i^^ . These fields describe 
the forces on the electrons arising from Berry phases 
which they pick up when their spin adiabatically follows 
Jl(r,i). They couple to the spin rather to the charge: 
electrons with magnetic moment parallel (antiparallel) 
to fl carry the "emergent charge" —1/2 (-1-1/2), respec- 
tively. For Vs ~ the change of the free energy density 
is given by 



dtF = ^ dtCi 
6n 



-aM{dt(lf - a'MiE^y 



(4) 



Mn X a' 



E^^+^{vsXB-)^^d,n. (3) 



which shows that the last term describes the dissi- 
pated power (X [E'^Y arising from the emergent electric 
field, a' M is therefore approximately given by the spin- 
conductivity CTs- 

We have also added a new /3'-term. The presence of 
such a term becomes evident if one considers the special 
ease of a Galilean invariant system. In this case, all forces 
have to cancel when the magnetic structure is comoving 
with the conduction electrons, d{r, t) = tl{r — Vst). This 
is only possible for a = /3 and a' = /3'. Solids are not 
Galilean invariant and therefore /3' is different from a' 
but one can, nevertheless, expect that the two quantities 
are of similar order of magnitude. 

Which of the d amping terms will dominate? As 
pointed out in Refs. (20ll2ll |. the naive argument, that the 
a' terms are suppressed compared to the a terms as they 
contain two more derivatives, is not correct. The dis- 
tance of skyrmions is^ proportional to 1/Aso, where Xso 
parametrizes the strength of spin-orbit coupling. While 
the a' term has two more gradients compared to the a 
term, the contribution arising from a' is, nevertheless, 
of the same order in powers of Xso, if "we assume that 
a arises only from spin-orbit coupling, a oc Agg, while 
a' (X AgQ (ohmic damping (see above) does not require 
spin-orbit effects). As furthermore a is proportional to a 
scattering rate while a' is proportional to a conductivity 
and therefore the scattering tim(M^, a' and /3' might 
be the dominating damping terms in good metals. 



III. DYNAMICS OF SKYRMIONS 

Our goal is to describe both the drift and the rotation 
of the skyrmion lattice in the limit of small current densi- 
ties and small magnetic or thermal gradients. We there- 
fore assume that is small compared to all characteristic 
velocity scales of the skyrmion lattice (e.g. Tc — T multi- 
plied with the skyrmion distance). The gradients should 
be so small that the total change across a domain of ra- 
dius rd remains small, VA <C A where A is -B or Tc — T 
for magnetic or thermal gradients, respectively. In this 
limit, both the drift velocities < Vs and the angular 
velocity dt4> oc ■ VA characterizing rotational motion 
remain small. Below we will show, that even dtiprd, the 
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velocity at the boundary of the domain remains smaU in 
the considered hmit. 

We can therefore neglect macroscopic deformations of 
the magnetic structure and consider the following ansatz 

n{r, t) = • Oo [R-^I) ■ (r - Vdt)) (5) 

Here flo{r) describes the static skyrmion lattice, is 
a matrix describing a rotation by the angle (jj around the 
direction of the skyrmion lines (i.e. around the field di- 
rection when anisotropics are neglected, which will be 
assumed in the following) and v^t describes the loca- 
tion of the center of the skyrmion domain. This ansatz 
describes a magnetic domain which rotates around its 
center, while the center is moving with the velocity Vd- 
When the torque forces are too weak to induce a steady- 
state rotation, such that dt(t> = 0, we will study rotations 
by the finite angle as in the experiment of Rcf . Q . 

A. Drift of domains 

To obtain an equation for the drift velocity Vd we fol- 
low Thielei^ and project Eq. ([3]) onto the translational 
mode by multiplying Eq. ([3]) with 9;$! and integrating 
over a unit cell (UC). We thereby obtain to order (VA)° 
(where no rotations occur) an equation for the force per 
2d magnetic unit cell (and per length)^ 

Gx{v.,- Vd) + T)0Vs - &Vd) + Fpin = (6) 
= / SrMBl^QB,, Q = AttMW 

JVC 

T>,j^ / (frMd^ndjCl^VPij 

JVC 

V = ( SrM{B^f 

JVC 

a = a + a'V' /V and /3 = ^ P'V' /V 

Here the first term describes the Magnus force which is 
proportional to the topological winding number W which 
is for the skyrmion lattice exactly given by W = — 1. 
G is called the gyromagnetic coupling vector following 
Thielei^. The second term are the dissipative forces with 
the projector P into the plane perpendicular to B, i-' = 

(l-B-B^). 

Besides the forces discussed above, also pinning forces, 
described by the last term in Eq. (jS]), have to be con- 
sidered. Formally, they are encoded in spatial fluctua- 
tions of SF/SCt in Eq. (jSj. The Thiele approach, used 
above, which considers only a global shift (or a global 
rotationi^, see below) of the magnetic structure does not 
capture these pinning effects as for a perfectly rigid mag- 
netic structure, random pinning forces average to zero, 
such that no net effect remains in the limit of a large 
domain. To describe pinning, it is necessar y to take 
into account that the magnetic structure adjusts locally 



to the pinning forces, a complicated problem for which 
presently no full solution existsiii^^ and which is far be- 
yond the scope of the present paper. Instead, we use a 
phenomenological ansatz and write for a finite drift ve- 
locity Vd 

Fpin = -47rAfWpin/(wd/Wpin) *d (7) 

to describe a net pinning force, which is oriented op- 
posite to the direction of motion. Its strength, which 
depends both on the number (and nature) of defects re- 
sponsible for pinnning and the elastic properties of the 
skyrmion lattice, is parametrized by the 'pinning veloc- 
ity' Wpin. The function /(x) with f{x — > 0) = 1 and 
f{x — !■ oo) = parametrizes the non-linear dependence 
of the pinning force on the velocity. Presently, it is not 
clear to what extent f{x) depends on microscopic details 
and also the exponent ly is not known. For large driving 
velocities, however, pinning becomes less and less impor- 
tant {v < I) i6ii7,29 ^ j£ ^jjg driving forces are smaller 
than the force AirMvpin, needed to depin the lattice, Vd 
vanishes and the pinning forces cancel exactly the driv- 
ing forces. Note that we do not consider creep, i.e. a 
tiny motion driven by thermal (or quantum) fluctuations, 
which occurs even in the pinning regime^. If the dissi- 
pative forces can be neglected, it is in principle possible 
to obtain f{x) from a measurement of the velocity of the 
skyrmion lattice^i. 

In the limit ^ 'fpin, where i^pin can be neglected, 
we solve Eq. ([6]) for Vs -L B to obtain 

^ + a3(p/g)^ + a r- + " g ^ " "-0 

with & = a + a'V /V and /3 = /3 + jS'V /V . 



B. Rotational torques 

By symmetry, a small uniform current cannot induce 
any rotational torques on a skyrmion lattice with per- 
fect sixfold rotation symmetry and therefore all effects 
arise from gradients. To derive an equation for the ro- 
tational torques which determine the rotations around 
the B axis, we foUowi^ a similar procedure as used for 
the translations by multiplying ([3]) by the generator of 
rotations applied to Ct 

d^n = B xn- {B{Ar xV))n (9) 

with Ar = r—Vdt and integrating over r. This procedure 
leads to several types of contributions. 

For the first type of contribution, we observe that the 
second term in Eq. ([3]), linear in Ar, is much larger than 
the first one which we can therefore neglect whenever 
the second term contributes. The second term induces 
torques of the form r x f where the force fi is obtained 
by multiplying Vift with the terms of Eq. ([3]). In the 



presence of gradients of the parameter A we obtain 

B-[rx /(A(r))] « J (^B ■ [r x dxf]) {r ■ VA) 

^A.B-[V\xdxjf] (10) 

where A is the area of the domain. Here it is essential to 
take the derivative with respect to A for fixed Vd reflect- 
ing that due to the rigidity of the skyrmion crystal Vd is 
constant across the domain. As the sum of all relevant 
forces vanishes [Eq. ©], J2i fii^^'^d) = 0, one obtains 
ilE,/, = while i:E^f^l, is finite. In Eq. m 
we have implicitely assumed a symmetrically shaped do- 
main, where integrals odd in r vanish. In general, there 
will also be a shape dependent torque T^hape arising even 
in the absence of a gradient. As its sign is random, it 
can easily be distinguished from the other torques (and 
appears to be relatively small in the MnSi experiments^) . 
More difhcult is the question what happens at the inter- 
face of different domains or when a domain comes close 
to the surface of the sample. Nominally surface forces are 
suppressed by a factor proportional to 1/ ^/A compared 
to the bulk terms considered above but the relevant prcf- 
actors are difficult to estimate. We will neglect in the 
following formulas both extra surface forces and shape 
dependent torques. 

A different contribution arises from the time deriva- 
tives dtfl = dt(l)d^Cl - {vdV)n in Eq. ©. The contri- 
bution proportional to Vd is of the form discussed above. 
The term proportional to dt4> leads to extra torques inde- 
pendent of VA. By combining the linear term in Ar from 
diptt with the second term of Eq. © we obtain for exam- 
ple the contribution adfcj) J M[{B[Ar x V])r2]^ which is 
also linear in A. Physically this term describes the fric- 
tional torque which is linear in the angular velocity dt(l). 
The frictional torque per volume is proportional to A be- 
cause the velocity and therefore the frictional forces grow 
linearly with the distance from the center of the rotating 
domain. 

Finally, a contribution exists which is independent of 
the gradients VA, the angular velocity 9(0 and of Vs- 
This contribution describes that in the absence of any ex- 
ternal perturbation the skyrmion lattice has a preferred 
orientation relative to the atomic lattice. Such terms ex- 
press that angular momentum can be transferred directly 
from the skyrmion lattice to the underlying atomic lat- 
tice mediated by spin-orbit coupling and small anisotropy 
terms (not included in Eg. ([T 1)V These terms have been 
discussed in detail in Ref. jl8| . This torque per unit cell 
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can be expressed by the change of free energy per unit 
cell, Ftjc, upon rotation by the angle </), where = 
reflects the equilibrium position and sin 60 reflects the 
sixfold symmetry of the skyrmion lattice. As has been 
discussed in Ref. [l^ , the absolute value of x in materials 



15 
10 

5 

^ 
-5 
-10 
-15 



1 1 


1 1 1 


1 1 1 


1 i_ 






A. 








— coy / 




" J 

V 

"I 1 


, 1 , 


, 1 , 


1 i" 



1.5 
1 

0.5 

,f 

-0.5 
-1 

-1.5 



-1.5 -1 -0.5 0.5 1 1.5 

VA,-V 



Figure 2; Rotation angle (in units of 1°) and angular ve- 
locity tj (times the prefactor 7) as a function of VA • de- 
termined from Eq. (|13p . 



like MnSi is tiny as it arises only to high order in spin- 
orbit coupling and, in contrast to all other terms, it is not 
linear in the size of the domain. Nevertheless, we have 
to consider this term, as it is the leading contribution 
arising to zeroth order in VA and Vg. 

Balancing all torques (per unit cell) we obtain as our 
central result 
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The direction of the torques, which depends on the rel- 
ative orientation of velocities and currents, is for Vd — 
(and dt(j) = 0) fully consistent with the simple picture 
shown in Fig. [TJ the dissipative torques To arise when 
gradient and current are perpendicular to each other 
while the reactive torque Tg arising from the Magnus 
force is activated for a parallel alignment of gradients 
and currents. For finite Vd, however, this simple intu- 
itive picture cannot be used especially as some of the 
torques tend to cancel when Vd approaches Vg- 



C. Rotation angle and angular velocity 



Eq. (|12p can be rewritten in the compact form 

sin60 = -7 9t0-h VA-V; (13) 

where 7 = and the vector Vs = Vs[vs] can be ob- 
tained by first solving Eq. ([5]) to obtain Vd as a function of 



Vs- This function is inserted into Eq. p2|) which, finally, 
is devided by — x. The function ys[i's] with Vs[0] ~ is 
proportional to the area A of the domain and encodes all 
information how the current couples to small gradients 
and includes contributions from Magnus forces, dissipa- 
tive forces and pinning. 



1. Dependence on size of gradients 

Qualitatively, three different regimes have to be dis- 
tinguished. For j < jc, when pinning forces cancel all 
reactive and dissipative forces, there is neither a motion 
nor a rotation of the skyrmion lattice, Vs = 0,0 = 0, 
within our approximation. Note, however, that it is 
well known from the physics of charge density waves or 
vorticesii that even below jc a slow creep motion is pos- 
sible. Whether during this creep also rotations are possi- 
ble is unclear, but the rather sharp onset of the rotation 
in the experiments of Ref. [13, see Fig. |8l seems to con- 
tradict a scenario of pronounced rotations during creep. 
For j > jc, the domains move and Vs will generally be 
finite. In this case, one can control the size and direction 
of rotations by the size of VA as shown in Fig. [21 For 
I VA -Vsl < 1, one obtains a solution where dtcp = but 
the gradients induce a rotation by a finite angle 



- arcsin VA • V, 
6 



(14) 



which grows upon increasing VA from zero until it 
reaches the maximal possible value 7r/12 = 15° (rotations 
by an average angle of 10° have already been observedi^, 
see Fig. [5]). For jVA-Vsl > 1 the domain rotates (see 
Fig. [2]) with the (average) angular velocity 
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Figure 3: Inset: Rotation angle (in units of 1°) as a function 
of time for three values of VA > 1, see Eq. (|16fl . For 
torques close to the value where rotations sets in, the rotation 
is strongly anharmonic. This can also be seen by considering 
the Fourier coefficients Cn = | J^^^'^'^ dtcpe'"^"* dt\ shown in 
the main panel as a function of VA Vs. 



For fixed uj the velocities at the boundary of the do- 
main, Vh = tdrd, grow linearly with the radius of the do- 
main rd- As we have assumed that the gradients across 
the sample and therefore also across a single domain are 
small, rd VA ^ A, the velocities nevertheless remain 
small, fb <C |Vs|A/7 < Vg/a. While our estimate does 
not rule out that Ub can become somewhat larger than Vs 
or Vd, we expect that the typical situation is that the ve- 
locity Wb arising from the rotation remains smaller than 
the overall drift velocity of the domain Vd- This estimate 
also implies that violent phenomena like the breakup of 
domains due to the rotation will probably not occur. 



and Eq. ([T^ is solved by 



4)(t) = — arctan 
3 



1 + juj tan(3wt) 



(16) 



displayed in the inset of Fig. |31 As both 7 and Vg are 
linear in the area A of the domain, w « (VA • Vs)/7 be- 
comes independent of the domain size for A 00. In this 
limit, the domain rotates continuously, (j) = tdt. Close to 
the threshold, VA • Vs = 1, however, the rotation be- 
comes very slow close to an angle of 15° (plus multiples 
of 60°). 

A way to detect the rotation of the magnetization is 
to exploit the emergent electric field E"^ which obtains 
a contribution proportional to oc dt(j) and can be mea- 
sured in a Hall experimental. In Fig. [3] we therefore 
show the modulus of the Fourier components, |c„| = 
I / e'6'^"*at</)dt| of dtcj) as a function of VA - Vs. At the 
threshold, all Fourier components are of equal weight 
while for large gradients the rotation gets more uniform. 



2. Domain size dependence and angular distribution 

In a real system, there will always be a distribution 
of domain sizes A. Both Vs and 7 are linear in A and 
therefore both the rotation angle p4)) and the angular 
velocity psp will in general depend on the domain size 
and therefore on the distribution of domains. 

Only in the limit |VA - Vsl ^ 1, the dependence on A 
cancels in Eq. ([T5)) and all domains rotate approximately 
with the same angular velocity. For | VA • Vs| ^ 1 one will 
in general obtain a distribution of rotation angles which 
can be calculated from the distribution of domain sizes 
Pd{A). For the static domains only angles up to 15° are 
possible with 



ps _ 



dAPd{A)S{^- 



arcsin(A/^c) 
6 



6Ac cos(60) Pd{Ac sm{6(j))) for < </> < — (17) 
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Figure 4: Angular distribution of the rotation angle of the 
skyrmion lattice for various values of Ao/Ac oc VA (see text). 
Here we assumed a distribution of domain sizes of the form 
Pd{A) = e"'*^^" While static domains contribute only for 

< < 15° one obtains a smooth angular distribution when 
one takes the rotating domains with A > Ac into account. 



where A^. = A/^VX • Vg) is the size of a 'criticar domain 
which just starts to rotate continuously. 

The continuously rotating domains also have a non- 
trivial angular distribution as their rotation will be 
slowed down when the counterforces are strongest, i.e., 
for cj) = 15°, see inset of Fig.[3l The angular distribution, 
P^, of the rotating domains is calculated from distribu- 
tion of domain sizes, Pd{A), and the angular distribution, 
p^{A)^ of a single domain 



pr _ 



dAP4A)p^^{A) 



Tdr 



3 ^A^-Al 
TT A — Ac sin 6(j) 



(18) 



where T = 27t/{6uj). While both P^ and P^ arc non- 
analytic at = 15°, the total distribution, P^ = P? + PI 



is smooth for > and normalized to 1, Jq"^^^ P<j, d(j) 



1. 



In Fig. U we show assuming the domain distribution 
PdiA) = e'^/^°4i for various values of AqMc- 

In elastic neutron scattering, the skyrmion phase is 
observed by six Bragg spots forming a regular hexagon in 
a plane perpendicular to the magnetic field. A rotation of 
the skyrmion domain results in a rotation of these Bragg 
spots. Therefore the angular distribution of rotation 
angles is directly observable (see Sec. IIIIDI below) by 
measuring the scattering intensity as a function of angle. 
By comparing angular distributions for different strength 
of the current or gradient, one can - at least in principle 
- obtain not only Ac as a function of VA or j but also the 



distribution of domain sizes. The latter can be extracted 
most easily in the regime where most of the domains 
do not rotate continuously by plotting P^/ cos 60 as a 
function of sin 60 using Eq. ([T7)) . 



3. Dependence on strength of current 

While the behavior of 4> and w as a function of VA is 
rather universal and independent of microscopic details, 
its dependence on the strength of the current for fixed 
VA is much more complex. As discussed above, Vg = 
for j < jc- Directly at jc, when the domain starts to 
move with ii^ ~ 0, Vs jumps to the finite value 



I Vs —V^ 




(19) 



Note that the jump is independent of a and a' as well as 
of their gradients, as the skyrmions are not moving di- 
rectly at the depinning transition (see Fig.|6]). Depending 
on the direction and size of VA, the jump of Vs either 
leads to a jump of the rotation angle for | VA • Vs| < 1 or 
immediately to a continuous rotation for | VA • l^j > 1. 

Upon increasing the current, VA • Vs can either in- 
crease, decrease or even change its sign depending on 
(i) the direction of VA and (ii) on the question which 
of the forces changes most strongly when varying A (i.e., 
temperature or magnetic field). 

Motivated by existing experimental data (discussed be- 
low in Sec, nil D|) we study the case of a temperature gra- 
dient, X = t, based on the following assumptions. First, 
we assume that all damping constants are temperature 
independent (this assumption is relaxed later). Second, 
we need also a theory for the temperature dependence of 
the pinning force. Here we use the experimental obser- 
vational that the critical current is almost temperature 
independent at least for a certain range of temperatures. 
Within our theory, Eqs. (|6]) and this implies that 
all temperature dependence of i^pin (i.e., the dependence 
on the parameter t in Eq. ([1])) arises from the tempera- 
ture dependence of the magnetization M which we cal- 
culate from the Ginzburg-Landau theory ([1]). From the 
Ginzburg-Landau theory, we obtain also the temperature 
dependence of the other parameters, see Fig. [T] 

In Fig. [5] we show a typical result (for temperature- 
independent dissipation constants) for the rotation angle 
and angular velocity of a skyrmion domain as a function 
of Vs in the presence of a temperature gradient. For a 
temperature gradient perpendicular to the current (lower 
panel of Fig. [S]), the rotation angle increases after the 
initial jump. For the gradient parallel to the current, 
however, we obtain that the rotation angle drops after 
the initial jump (upper panel). For larger values of Vg 
the angle rises again until it reaches its maximal value of 
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Figure 5: Rotation angle cj) (in units of 1°) and angular veloc- 
ity, 'y(jj, as a function of Vs for a temperature gradient parallel 
(Vt||vs, — (—0.1,0,0), upper panel) and perpendicular 
(Vt _L Va, Vf — (0,-0.05,0), lower panel) to the current 
(a = 0.2, /3 = 0.45, a' = 0.01, /3' = 0.2, A/x = 200, t = 
-1,B = (0,0, l/^/2),l;pin = 1,/ = 1). For both geometries 
one observes a jump of at from zero to a finite 

rotation angle. After the initial jump the rotation angle in- 
creases for the perpendicular configuration (panel b) while for 
the parallel arrangement first a drop and then an increase up 
to the maximal angle of 15° occurs. For larger Vs a contin- 
uous rotation characterized by the angular velocity lj sets in 
for both configurations. For the calculation we assumed that 
the damping parameters and Vpin are independent of t. 



15°. This qualitative shape of the curve appears to be 
rather independent of the precise values of the various 
parameters if we assume that all damping parameters 
are temperature independent. 

In Fig. [Sj we plot the rotation angle for small cur- 
rent densities taking an extra effect into account which is 
present in the experiments described in Ref. [l^: as the 
temperature gradients are induced by the currents, they 
grow quadratically with Vs . This does not give rise to any 
qualitative changes. The thin blue curve Fig. [Sj thereby 
reflects the same physics as the corresponding curve in 



Figure 6: Rotation angle (fi (in units of 1°) as a function of Vs 
for a temperature gradient parallel (Vt||vs, upper panel) and 
perpendicular (Vt _L Vs, lower panel) to the current. The pa- 
rameters are the same as in Fig. [S] with two exceptions. First, 
we have taken into account that in the experiments of Ref. jTy 
the temperature gradient grow with the square of the applied 
current, Vt = (-0.1w?,0,0) and Vt = (0, -0.051;?, 0)), for 
current parallel an perpendicular to Vs, respectively. For the 
thin blue curve we assumed (as in Fig. that the damping 
constants are independent of t while for the thick green curve 
a weak temperature dependence of the damping constant a, 
Va = 0.035 Vt, was assumed. This parameter has been cho- 
sen to reflect the experimental observation, see Fig. [S] For 
even stronger currents (not measured experimentally and not 
shown in the figure) the size of the torque drops again and a 
finite rotation angle is obtained for 1.57 ^ Vs ^ 2.53 in the 
parallel configuration with the temperature dependent damp- 
ing constant. 



Fig. E] (note the different scale on the x axis) . The thick 
green curve of Fig. |6| shows that one can, however, ob- 
tain qualitative different results (an increase rather than 
a reduction of the rotation angle after the initial jump 
for T gradients parallel to the current, upper panel) by 
including a small temperature dependence of the Gilbert 
damping a. As we will discuss in Sec. IIIl Dl this can 
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in which the skyrmion lattice drifts. 

Vj {v,~Vd)\\ 



(20) 



The ratio can be obtained experimentally by mea- 
suring the rotation angle or the angular velocity for VA 
parallel and perpendicular to the current, from which 

one can obtain directly using Eqs. and (|15p . 
For small angles, arcsina; « x, for example, one obtains 
directly from the ratio of the two rotation angles. A 
different, but probably more precise way to determine 
this ratio is to find experimentally the "magic angle" 
(pm of gradient vs. current, where all rotations vanish, 
VA • V, = 0. In this case one obtains 



Figure 7: Change of C/,© and V defined in Eq. (O with 
temperature t. The applied mag netic field is /i/V2(0,0, 1). 
Dashed lines are for h — 0.9 and continuous lines for h = 1.1. 



reproduce qualitatively the experimentally observed be- 
havior. 



4- Dependence on orientation of gradients 

Fig. [5] shows that the rotational torques on the sys- 
tem depend strongly on the relative orientation of gradi- 
ent and current. More importantly, one probes different 
physical mechanism for gradients parallel or perpendicu- 
lar to the current. This effect was already discussed in the 
introduction, see Fig. [TJ where, however, only the simple 
case of a static domain without pinning was described. 
In reality, the situation is more complex. All directional 
information is encoded in the function Vs{vs) which can 
be obtained by first solving Eq. (jHl) to obtain Vd and then 
comparing Eqs. and ([T^ . Unfortunately, a rather 
large number of unknown parameters (most importantly, 
the pinning forces and their dependence on A) enters the 
description. Therefore we will discuss in the following 
only a few limiting cases. 

A drastically simplified picture occurs in regimes when 
only two forces dominate in Eq. ([S]). For example, close 
to the pinning transition, the Magnus force is of the 
same order as the pinning force while the two dissipa- 
tive forces are typically much smaller. In this case one 
can use Eq. ^ to show that Vd becomes proportional 
to B X {vs — Vd). Thus, for an A-independent Vg, both 
the reactive rotational coupling vector and the rotational 
pinning vector become proportional to VA • {v^ — Vd) 
(here we neglect a possible A-dependence of v^). There- 
fore the ratio of the component of Vs parallel (Vs") and 
perpendicular (Vs"*") to Vg depends only on the direction 



Vs- 



1 



tan( 



(21) 



This should allow for a quantitative determination of 
/ ujj . As v\ can be measured independently us- 
ing emergent electric fields generated by the motion of 
skyrmions^l, one can obtain the complete information 
on the drift motion by combining both experiments. It 
is also instructive to compare skyrmions and vortices in 
a superconductor. Vortices and skyrmions follow essen- 
tially the same equation of motions, Eq. The rele- 
vant parameters (and therefore also the pinning physics) 
are, however, rather different. For vortices in conven- 
tional superconductorsii^ the dissipation is very large 
"Da ^ Q. Therefore, vortices drift - up to small correc- 
tions - predominantly perpendicular to the current while 
for magnetic skyrmions we expect that at least not too 
close to the depinning transition, the motion is domi- 
nantly parallel to the current. 

In the limit where the pinning forces can be neglected, 
i.e., Vs 3> Wpin; to linear order in /3 and a the vector Vg 
is given by 



A 
Attx 



A 
Attx 



B X V 



B X V. 




dV0 




Here we also neglected a possible A-dependcncc of Vs ■ In 
this limit the rotation can be induced primarily by gra- 
dients perpendicular to Vs reflecting that the motion of 
skyrmions is mainly parallel to the current, see Eq. ((20)) 
and Eq. ([8|). This is also consistent with the behavior 
shown in Fig. [5] where we used a two-times smaller gradi- 
ent for the perpendicular configuration and obtained nev- 
ertheless an onset of the rotational motion for values of 
Vs much smaller than in the parallel configuration. Note 
that in a Galilean invariant system, a — (3, no torques 
can be expected. 
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D. Experimental situation 

Our study is directly motivated by recent neutron 
scattering experiments in the skyrmion lattice phase of 
MnSii^. In the presence of a sufficiently large current, 
a rotation of the magnetic diffraction pattern by a finite 
angle was observed when simultaneously a temperature 
gradient was present (only temperature gradients paral- 
lel to the current have been studied). The rotation angle 
could be reversed by reversing either the direction of the 
current, the direction of the magnetic field or the direc- 
tion of the temperature gradient. This clearly showed 
that rotational torques in the experiment were driven by 
the interplay of gradients and currents as studied in this 
paper. 

In Fig. [HK we reproduce Fig. 3 (A) of Rcf. [l3l, which 
shows the average rotation angle (defined as the maxi- 
mum of the azimuthal distribution of the scattering in- 
tensity) as a function of current density. Above a critical 
current, j > jc, the rotation sets in. The rotation angle 
initially increases abruptly, followed by a slower increase 
for larger current densities. When comparing these re- 
sults with our theory one has to take into account that 
the temperature gradient in the experiment was not inde- 
pendent of the strength of the applied electrical current 
density as it originated in the resistive heating in the 
sample. Therefore the temperature gradient was grow- 
ing with (i.e. the heating rate due to the electric cur- 
rent). This was taken into account in Fig. [H] as discussed 
above. For a full quantitative comparison of theory and 
experiment, it would be desirable to have data, where 
the applied current as well as both the strength and the 
direction of the gradients are changed independently. As 
such data is presently not available, we restrict ourselves 
to a few more qualitative observations. 

In our theory wc expect a jump of the rotation angle 
at jc, which depends on the domain size. This appears 
to be consistent with the steep increase of the rotation 
angle as observed experimentally at jc, especially when 
taking into account the experimental results are subject 
to a distribution of domain sizes. 

Interestingly, the experimentally observed increase of 
the rotation angle after its initial jump is apparently 
not consistent with the predictions from the extended 
Landau-Lifshitz-Gilbert equation shown in Eq. ^ if we 
assume a, a', /?,/?' are independent of temperature. As 
shown in Fig. [6l we can, however, describe the experi- 
mentally observed behavior if we assume a weak temper- 
ature dependence of the Gilbert damping. 

An important question concerns, whether the existing 
experiments already include evidence of some larger do- 
mains that rotate continuously. Fig. [5^ shows that for the 
largest currents average rotation angles of up to 10° have 
been obtained. As this is rather close to the maximally 
possible value of 15° for static domains, this suggests that 
continuously rotating domains are either already present 
in the system or may be reached by using slightly larger 
currents or temperature gradients. 



We have therefore investigated the angular distribu- 
tion of the scattering pattern usingthe same set of ex- 
perimental data analyzed in Ref. [ifll (technical details 
of the experimental setup are reported in this paper) . In 
Fig-Eb we show the azimuthal intensity distribution with 
and without applied current. Already for zero current 
a substantial broadening of the intensity distribution is 
observed. The origin of this broadening are demagneti- 
zation effects which lead to small variations of the orien- 
tation of the local magnetic fields in the sample tracked 
closely by the skyrmions. It has been shownii that this 
effect can be avoided in thin samples when illuminating 
only the central part of this sample. For the existing 
data this implies that a quantitative analysis of is not 
possible. We observed that the measured experimental 
distribution of angles extents up just to 15°. Therefore, 
from the present data we can neither claim nor exclude 
that continuously rotating domains already exist for this 
set of data but slightly larger current densities or gradi- 
ents should be sufficient to create those. 



IV. CONCLUSIONS 

The magnetic skyrmion lattices, first observed in MnSi, 
have by now been observed in a wide range of cubic, 
chiral materials including insulators^Sii^, doped semicon- 
ductors^ and good metals^ii^^. This is expected from the- 
ory: in any material with B20 symmetry, which would 
be ferromagnetic in the absence of spin-orbit coupling, 
weak Dzyaloshinskii Moriya interaction induce skyrmion 
lattices in a small magnetic field. While in bulk they are 
only stabilized in a small temperature window by ther- 
mal fluctuations close to the critical temperature, they 
are much more stable in thin films^i^^. 

From the viewpoint of spintronics, such skyrmions 
are ideal model systems to investigate the coupling of 
electric-, thermal- or spin currents to magnetic textures: 
(i) the coupling by Berry phases to the quantized wind- 
ing number provides a universal mechanism to create effi- 
ciently Magnus forces, (ii) skyrmion lattice can be manip- 
ulated by extremely small forces induced by ultrasmall 
currentsi2i21, (iii) the small currents imply that also new 
types of experiments (e.g., neutron scattering on bulk 
samples) are possible. 

We think that the investigation of the rotational dy- 
namics of skyrmion domains provides a very useful 
method to learn in more detail which forces affect the dy- 
namics of the magnetic texture. As we have shown, the 
rotational torques can be controlled by both the strength 
and the direction of field- or temperature gradients in 
combination with electric currents. They react very sen- 
sitively not only on the relative strength of the various 
forces but also on how the forces depend on temperature 
and field. 

While some aspects of the theory, e.g. the dependence 
on the strength of the gradients, can be worked out in 
detail, many other questions remain open. An important 
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question is, for example, to identify the leading damping 
mechanisms and their dependence on temperature and 
field. Also an understanding of the interplay of pinning 
physics, damping and the motion of magnetic textures 
is required to control spin torque effects. Here future 
rotation experiments are expected to give valuable in- 
formation. Furthermore, it will be interesting to study 
the pinning physics in detail and to learn to what ex- 
tent skyrmions and vortices in superconductors behave 
differently. 

One way to observe the rotation of the skyrmion lattice 
is to investigate the angular distribution of the neutron 
scattering pattern as discussed in Sec. IIII CI This does, 
however, only provide indirect evidence on the expected 
continuous rotation of the skyrmion lattice. Therefore 
it would be interesting to observe the continuous rota- 
tion more directly. For example, one can use that time- 
dependent Berry phases arising from moving skyrmions 
induce "emergent" electrical fields which can be directly 
measured^l in a Hall experiment. Here it would be in- 
teresting to observe higher harmonics in the signal which 
are expected to appear close to the threshold where con- 
tinuous rotations set in, see Fig. |31 

In future, it might also be interesting to use instead of 
electrical current other methods, e.g. pure spin currents 
or thermal currents, to manipulate skyrmion lattices (e.g. 
in insulators). We expect that also in such systems the 
investigation of rotational motion driven by gradients will 
give useful insight in the control of magnetism beyond 
thermal equilibrium. 



Figure 8: a) Average rotation angle Acfi (in units of 1°) of the 
skyrmion lattice in MnSi measured by neutron scattering in 
the presence of an electric current and a temperature gradient 
parallel to the current. The figure is taken from Ref. [l^] 
where further details on the experimental setup can be found, 
b) Angular distribution of the intensity normalized to 1 
for currents of strength j = (black diamonds) and j ~ 
-2.07- 10*^ A/m^ for T = 27.4 K (red circles). The lines are 
Gaussian fits serving as a guide to the eye. The distribution of 
angles extents up the maximally possible rotation angle of 15° 
which suggests that some of the larger domains are rotating 
with finite angular velocity for this parameter range. 
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